Mixing Times in Quantum Walks on the Hypercube 
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The mixing time of a discrete-time quantum walk on the hypercube is considered. The mean 
probability distribution of a Markov chain on a hypercube is known to mix to a uniform distribution 
in time 0(n log n). We show that the mean probability distribution of a discrete-time quantum 
walk on a hypercube mixes to a (generally non-uniform) distribution ti(x) in time 0(n) and the 
stationary distribution is determined by the initial state of the walk. An explicit expression for 
■k{x) is derived for the particular case of a symmetric walk. These results are consistent with 
those obtained previously for a continuous-time quantum walk. The effect of decoherence due to 
randomly breaking links between connected sites in the hypercube is also considered. We find that 
the probability distribution mixes to the uniform distribution as expected. However, the mixing 
time has a minimum at a critical decoherence rate p ~ 0.1. A similar effect was previously reported 
for the QW on the A^-cycle with decoherence from repeated measurements of position. A controlled 
amount of decoherence helps to obtain — and preserve — a uniform distribution over the 2" sites of 
the hypercube in the shortest possible time. 

PACS numbers: 03.67.Lx, 05.40.Fb, 03.65.Yz 



I. INTRODUCTION 



Quantum walks (QW) are quantum versions of Markov 
chains in which coherent superpositions play a key 
role |1|. Both the continuous-time (CTC^W) and the 
discrete-time (DTQW) versions of QWs [2, 3] have been 
used as a basis for quantum algorithms that outperform 
their classical counterparts 0, H, 0, S HI • Some of them 
are based on a QW on the hypercube, which has received 
considerable attention after it was shown that the mean 
time required to traverse the hypercube between two op- 
posite vertexes (hitting time) is exponentially faster than 
its classical counterpart d, . In this work we consider 
the mixing properties of a DTQW in the hypercube, both 
in the coherent case and in the presence of decoherence 
from randomly broken links. 

The probability distribution of a classical Markov 
chain in a hypercube mixes to the uniform distribution 
in time 0(n log n), where n is the dimension of the hy- 
percube. In the quantum case, the probability distri- 
bution results from a unitary process and so it cannot 
converge to a stationary distribution. Aharonov et al 
have shown that the time-averaged distribution, defined 
in the discrete-time case by P{x,T) = y X^tllo^ ^(-^i 
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does converge to a limiting distribution One may 

sample from this distribution by selecting t uniformly in 
[0,T — 1], running the QW for t steps and then measur- 
ing the position of the walker. The time it takes for the 
average distribution to be e-close to this limiting distri- 
bution is known as the mixing time Me (we shall later 
provide formal definitions). Aharonov et al give an up- 
per bound for the mixing time of QWs on generic graphs. 
For a hypercube driven by a Grover coin, this bound is 
0(2— ). Our numerical simulations indicate a mixing 
time of 0{n/e), which is better than the classical mixing 
time. Moore and Russell have considered mixing times 
with respect to a uniform distribution in both the DTQW 
and CTQW on hypercubes 0. In the case of the CTQW 
they found that the average limiting distribution is not 
uniform and that one can always find e > such that 
no mixing time to the uniform distribution exists. Since 
the CTQW can be obtained from a DTQW by a suit- 
able limiting process we could expect the limiting 
distribution in the DTQW to be also non-uniform. 

In this work, we show that the asymptotic distribution 
of a DTQW on the hypercube is not uniform, obtain 
an explicit expression for it and characterize the average 
mixing time to this distribution. The instantaneous mix- 
ing time is a useful notion that captures the first instant 
in which the position distribution is e— close to a refer- 
ence distribution. It has been shown that, for both kinds 
of QW on a hypercube, the instantaneous mixing time to 
the uniform distribution depends linearly on the dimen- 
sion n [12]. This represents an improvement over the cor- 
responding classical walk, from 0(n log n) to 0{n). We 
confirm this result for the DTQW and show that e > 
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can be found such that an instantaneous mixing time to 
the stationary (non- uniform) distribution does not exist. 

All these properties are affected by decoherence. The 
impact of decoherence on QWs has been investigated 
mostly in one or two-dimensional systems using repeated 
measurements Il5| or topological noise from broken 
links Kendon and Tregenna have presented 

early numerical results showing the impact of repeated 
measurements on the hitting time (starting at node x, it 
is the mean time required to reach the opposite node, x) 
of a DTQW in a hypercube [l5S| . In the presence of deco- 
herence the natural expectation is that the position dis- 
tribution of a QW in an n-hypercube mixes to a uniform 
distribution essentially as in the classical case. Recently, 
the effect of repeated measurements on a QW on the hy- 
percube has been considered analytically by Alagic and 
Russell using superoperator techniques for the case of 
the CTQW [i^. They found that for weak decoherence 
rates, both the hitting time and the mixing time (to the 
uniform distribution) remain linear in n but if the rate of 
measurements is beyond a model-dependent threshold, 
the classical result is re-obtained, i.e. the CTQW dis- 
tribution mixes to uniform in time 0{nlogn). We con- 
sider here the effect on the probability distribution of the 
DTQW of randomly breaking links in the network. This 
mechanism involves no measurements and is an example 
of unitary noise [20| . 

Decoherence is usually regarded as an obstacle to quan- 
tum computing as it destroys the entanglement and su- 
perpositions required by a quantum processor. However, 
it has been argued that a controlled amount of deco- 
herence might be useful to obtain particular probability 
distributions. For instance, weak decoherence may be 
used to generate a quasi-uniform distribution in the QW 
on a line Maloyer and Tregenna have shown that 

the mixing time of a DTQW on an TV-cycle may be re- 
duced by allowing for some decoherence from repeated 
measurements, provided that those measurements affect 
the position of the QW ^21,]. We find a similar effect in 
the case of the hypercube with decoherence from bro- 
ken links. There is a critical decoherence rate for which 
the mixing time has a minimum. This result provides a 
new example in which decoherence, from topological ran- 
domness and in the absence of measurements, enhances 
a useful quantum feature, namely fast mixing times. 

This paper is organized as follows. In Section [ll] we 
consider the coherent QW on the hypercube and derive 
its asymptotic distribution which is non-uniform. In Sec- 
tion lllll we give the unitary evolution rule in the presence 
of broken links and present numerical results for the mix- 
ing time as a function of the decoherence rate. In Section 
IIVI we summarize our main results and present our con- 
clusions. 



II. COHERENT QW ON THE HYPERCUBE 

A coined quantum walk in a n-dimensional hyper- 
cube has a Hilbert space He Hp, where Tic is the 
n-dimensional coin subspace and Tip the 2"-dimensional 
position subspace. A basis for He is the set for 
< J < ?^ — 1 and Hp is spanned by the set {\x)} with, 
< a; < 2" — 1. A generic state of the discrete quantum 
walker in the hypercube is 

n-l 2"-l 

i*w) = E E ^Mj^^), (1) 

where \j,x) = \j) (8) |a;). The evolution operator for one 
step of the walk is 

U = So{C<»I), (2) 

where C = X)i j '-^ijN)OI is a unitary coin operation in 
He, I is the identity in Hp and S is the shift operator 
defined by 

n 2"-l 

^ = E E (3) 

j=0 x=0 

Here, 3C ® Cj is the bitwise binary sum between the n- 
component binary vectors, x = {xn-i, . . .xi, xq) and ej, 
a null vector except for a single 1 entry in the j*'' com- 
ponent. We will use both representations (decimal or 
binary) as needed. They may easily be recognized by the 
context in which they are placed. 

Applying the evolution operator on state ^ we obtain 
the map for the components of the wavevector, 

n-l 

ijt.xit + 1) = E V'j.^®e.(t)- (4) 

3=0 

The analysis of the problem is simplified in Fourier 
space. Since the hypercube is a Cayley graph of Zj, 
the most suitable transform is the Fourier transform on 
Z2 , spanned by the 2" kets 

|fc)^-l= ^(~i)^--|x), fce[o,2"-i], 

with bitwise product k ■ x = Y^"^q xjkj. The trans- 
formed amplitudes are V'l./c = Z]^=o^(~l)''''^V'i,a:- 
The evolution operator is diagonal in fc-space and acts 
non-trivially in the coin subspace. The operator Uk, 
which acts on = (^|^) = J27=o "^^M^) ^ matrix 
elements given by Uk{i,j) ~ (— 1) 'Cy. 

From this point on, let us particularize to the n- 
dimensional Grover coin, C{i,j) = 2/n — Sij, which 
obeys the permutation symmetry of the hypercube and 
is, among the operators of this kind, the one farthest 
from the identity We shall now describe the eigen- 
problem of Uk. Its eigenvalues depend only on n and 
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on the Hamming weight of k, defined as = jyj=o 
However, its eigenvectors, \vi{k)) for i — 1 . . . n, depend 
on k. 

For Hamming weights |fc| =0 (or = n), a set of 
71—1 degenerate eigenvectors with eigenvalue A = 1 (or 
-1) is |z/,(0)) = (|0) - for i e [l,n- 1]. The 

remaining eigenvector, with eigenvalue A = — 1 (or 1) is 
the uniform superposition |^'„(0)) = X]j=o^ b)- 

In the case where the Hamming weight takes values 
< |fc| < n, there are n — — 1 degenerate eigenvec- 
tors with eigenvalue — 1 , given by \vi{k)) = (|0) — \i))/y/2 
with i g [l,n— \k\ — 1] and |fc| — 1 degenerate eigenvec- 
tors with eigenvalue 1, given by \i'i{k)) — {\n — \k\) — 
|i 1))/V2 for i e [n - |fc|,n - 2]. The two remaining 
eigenvalues turn out to be the most relevant. They can 
be expressed as e'^^'^'' , with uj^ defined by 

coswfc = 1 . (5) 

n 

The corresponding conjugate eigenvectors are |^'„(fc)) = 

E"=o«i(^)li> and |i^„_i(fc)) = E"=oai(fc)lj) with 
components aj{k) given by 

In the last equation i = \/— T. This set of normal- 
ized eigenvectors, summarized in Table 1, forms a non- 
orthogonal basis. 

A. Limiting Distribution 

Let P{x, t) be the probability to find the walker on the 
vertex x of the hypercube at time t. As mentioned in the 
introduction, this probability depends on the initial con- 
dition and, as is typical of unitary evolutions, it does not 
converge to a stationary distribution. However, the time- 
averaged distribution P(x, T) = ^Y^J^q P{x,t), always 
converges as T — > oo [ll| . We define 

nix) ^ \im P{x,T) (7) 

i — >oo 

as the stationary distribution, and consider the problem 
of determining it for a QW in a hypercube of dimension 
n. The initial statevector is localized at vertex a; = and 
uniformly distributed in the coin subspace, 

n 

This choice respects the permutation symmetry of the 
hypercube. The initial state is expressed in terms of the 
eigenvectors of Uk as 

2"-l n 

|*(0)) = ^=^^a,(fc)|^,(fc))®|fc). 



where the coefficients are sums of components of the 
eigenvectors, ai{k) = ■^y=^ J2^Iq {i^i{k)\j) . These sums 
are zero unless i — n — 1 or i — n, so only these two 
eigenstates contribute, 

2"-l 

|*(0))= (a„_i(fc)|i.„_i(fc))+a„(fc)|i.„(fc)))®|fc). 

k=Q 

The relevant coefficients are 

a„(fc) = {y/\F\+iy/n- |fc|)/\/n2"+i and 
an-i(fc) = a^(fc) when |fc|G[l,n— 1]. If |fc| = or 
|fc| = n, then a„ = l/\/2"+i and a„_i — 0. The state 
of the walker at time t is 

2"-l 

\^{t))^ J2 (a„-i(fc)e-'"->„_i(A;)) + 

k=0 

an{k)e'^'''\iyJk)))<E>\k). (9) 

The probability of finding the walker at time t at vertex 
X is P{x,t) = J2]jZo IO'^l^(^))l^- This quantity can be 
evaluated in the fc-representation, 

2"-l n-l 

P(^'i) = F ^ (-l)('=®'=')-5](j-,fc|vl/(i))(*(i)|fc',j) 

k,k'=a j=0 
2"-l 

= — "Y (^_-[^-j(fcefc )-a; ^ 
^ k,k'=a 

{a„_i(fc)a;_i(fc')(i^„(fc)|j/„(fc')>e"'^"'=""^''* + 

a„(fc)a;(fc')(^^„-i(fc)|^^„-i(fc'))e'("'"-"'='^* + 
a„ (fc)a;_ 1 (fc') ('^n- 1 (fc) I '^n (fc')> e'^"* + 
an-iik)alik')(Mk)Wn-i{k'))e-'^^^+^'''^*}. 

Our goal is to calculate the asymptotic probabil- 
ity distribution tt{x) defined by Eq. ([7]). The 
first two terms contribute only if |fc| — \k'\, since 

liniT^oo T L^^o^ e"^*^"" = ^\k\,\k'\ and there is 
no contribution from the last two terms because 

limT^oo ELo^ e=^'^"''"^"'='-" = 0- Thus, the asymptotic 
distribution can be evaluated from 

k,k'=a 

{\an-i{k)\^Mk)Wn{k')) + 

|a„(fc)|'(^^„-i(fc)|;^„-i(fc')>}- (10) 

The overlap between the non-trivial eigenvectors is 

{Mk)Wn{k')) - {l^n-l{k)Wn-l{k')) 

^ n{k-k') + \k\{n~2\k\) 

2\k\{n-\k\) ■ ^ ^ 
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Table I: Eigenvalues and eigenstates of Uk- The quantities Uk and ctjik) are defined in Eqs. ([Sjl and respectively. 
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-f 
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After using |a„(fc)|2 = |a„_i(fc)|2 = 1/2'^ 
the stationary distribution 



we obtain 



2 
22^ 



1 

22^ 



2"-l 



E (-1) 



(|fe|=|fcVo,«) 



n(fc-fc') + |fc|(n-2|fe|) 
m{n~\k\) 



(12) 



The above expression is not efficient to calculate 'n{x). 
It is possible to simplify it by noting that ■k{x) depends 
only on \x\ = X]j=o After some algebra we eventually 
obtain 



f, -, n-1 i M W\ m 

1=1 j=0 /=0 m=Op=0 

|x| ^ m \ / m\ / n — \x\ — i + m 
J~ni + pJ\pJ\i~j — l + m — py 

— m\ /n — \x\\ l'\x\\ i{n~2i) + nj 

— p J \ i ~ m J \m J 2i {n — i) 



where the combinatorial coefficients are (^) 



(13) 



(n— m)! m! 

for n > TO > and (J^J = otherwise. This expression is 
equivalent to Eq. (jl2|) and, for some values of x, it yields 
simple results, such as 



1 



7^ (0) = 

\ ' An 



(14) 



It should also be noted that the identity 
7t{x) ~ 7r(2" — 1 — x) also helps in the evaluation 
of tt(x). Clearly, the asymptotic distribution for P[x,t) 
is not uniform for the initial condition of Eq.®. Note, 
for example, that for n sufficiently large, Eq. (fH)) leads 
to 7r(0) K. l/^27r(2n + 1) > 2"". We have performed 
numerical implementations which confirm the prediction 
of Eq. p^ . This distribution is shown for hypercubes 
of several dimensions in Fig. [TJ The maximum of the 
distribution occurs at the initial site Xo = and at 
xo = 2" — 1. Notice that ■k{x) takes only 1 + [n/2j 
different values which correspond to sites with Hamming 



distances 0, 1, 2 . . . [n/2j either to the initial site xq or 
to its opposite site, xq. 

Further insight may be obtained if one considers the 
probability p(|a;|) of finding the walker at long times at 
a Hamming distance from the initial site. There are 
(|"|) sites with Hamming weight |a;| and this probability 
is simply 



p{\x\) 



t:{x). 



(15) 



The stationary distribution ■k{x) depends only on \x\ and 
may be evaluated efhciently from Eq. lfT5|) . The probabil- 
ity p(|a::|) takes at most 1+ [n/2j different values, as 7r(x) 
does. In Fig. [2] J3(|a;|) is compared with the correspond- 
ing probability for a uniform distribution over the hyper- 
cube, in which case the most likely situation would be 
to find the walker with a Hamming distance ~ n/2 from 
the starting site. Instead, we see that the most probable 
situation is to find it near the marked sites xo , a;o , with 
\x\ ~ or |a;| n. 

In summary, the fact that the QW starts localized at 
xq = marks both xq and xq as special sites and per- 
sistent interference effects give rise to a non-uniform sta- 
tionary distribution over sites. The probability of find- 
ing the particle with a given Hamming distance varies 
much more slowly than for a uniform distribution over 
sites. We should remark that the results obtained for the 
asymptotic distribution are dependent on the initial con- 
dition. For instance, since the uniform superposition of 
coin and position eigenstates, -y=^ ^]=o Y^t:=o^ b' ^) ' 
an eigenstate of U with eigenvalue 1, for this particular 
initial condition the distribution remains uniform. 



B. Mixing time of a coherent evolution 

In this section we consider the mixing time for a coher- 
ent evolution. The rate at which the average probability 
distribution of a QW approaches its asymptotic distribu- 
tion is captured by the following definition [ll|. 

Definition II. 1 The average mixing time M,, of a quan- 
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Figure 1: Limiting distributions for QW in hypercubes with n = 3, 4, 6, 8 obtained from Eq. (|13p with the initial condition ([8}. 
As a reference we show the uniform distribution as a dashed horizontal line. 
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Figure 2: Asymptotic probability of finding the walker with a 
Hamming distance to the initial site, given by Eq. (|15p . for 
n = 25. The binomial distribution 2^(|"|), which corresponds 
to a uniform position distribution of the walker, is shown for 
comparison. 

turn Markov chain to a reference distribution tt is 

M, = min{T I Vi > T, \\Pt - 7r|| < e}, 

where \\A — B\\ = \ A(x) — B(x)\ is the total variation 
distance between the two distributions. 

An alternative definition captures the first instant in 
which the walk is e-close to the reference distribution 

Definition II. 2 The instantaneous mixing time of a 
quantum Markov chain is 

/, =min{i| ||Pt-^|| <e}. 



Both mixing times depend on the initial condition of the 
QW. 

The left panel of Fig. [3] shows the time dependence 
of the variation distance of the average distribution to 
the stationary distribution, Eq. (fT2|) . For long times, the 
variation distance decays approximately as ~ l/t while 
the corresponding distance to the uniform distribution 
remains essentially constant. 

For a quantum walk in a generic graph with arbitrary 
initial condition, an upper bound for the total variation 
distance to the asyrnptotic distribution Tr{x) was derived 
by Aharonov et al. 

||P(:z:,r)-7r(:z:)|i<^(l + ln(n2"-i)). 

This bound is expressed in terms of the minimum sepa- 
ration A between distinct eigenvalues of U . For a hyper- 
cube of dimension n driven by a Grover coin, this min- 
imum separation is A = min je*'^'' — 1| = 2/ ^/n. Thus, 
for large n, the average mixing time Mg is bounded by 

O {j^^T~^' right panel of Fig. [3] shows the linear 

dependence of the average mixing time with the dimen- 
sion n and the threshold e, for the initial state given by 
Eq. (IH). For sufficiently large dimension, it scales as n/e, 
which is consistent with this bound. 

The instantaneous mixing time, /g, the first time for 
which the position distribution P{x^ t) is e-close to a given 
distribution, must be calculated considering the appro- 
priate parity of each step of the walk. In the left panel of 
Fig. [4] we show the total variation distance of the instan- 
taneous distribution to the stationary distribution t:{x), 
given by Eq. (fT5|) . as a function of t/n, for several values 
of n. Note that the local minima do not approach zero 
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Figure 3: Left: total variation distance from the average distribution at time t to both the uniform and the stationary 
distributions, for a coherent QW moving on a hypercube of dimension n = 8. Axes are logarithmic in the larger plot and linear 
in the inset. Right: Mixing time as a function of the dimension n for several thresholds e. 
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Figure 5: Instantaneous mixing time to the uniform distribu- 
tion as a function of the dimension n. 



as t increases. Hence, for e < 0.3 and fixed n, there is no 
e-instantaneous mixing time. In the right panel of Fig. [4] 
we show the instantaneous mixing time as a function of 
the dimension n of the hypercube. There are large os- 



cillations in /e for small e and the dependence on n is 
clearly nonlinear. 

Moore and Russell explore the instantaneous mix- 
ing time to the uniform distribution with appropriate 
parity. This is a useful notion that captures the first in- 
stant in which the distribution of the QW is e-close to 
uniform. In Fig. [5] we show this e-instantaneous mixing 
time as a function of the dimension n of the hypercube, 
for several thresholds. It scales linearly with n with a 
slope close to tt/A. As reported in — and confirmed 
by our numerical calculations — , for t/n = 7r/4, the posi- 
tion distribution of the QW on the n-hypercube is close 
to uniform. However, the variation distance to the sta- 
tionary distribution ^{x), does not have this property, as 
shown in the left panel of Fig. 21 



III. DECOHERENCE AND MIXING TIMES 

Let us now consider the effects of decoherence on the 
stationary distribution and on the average mixing time 
of the QW on a hypercube. As a source of decoherence 
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we consider topological noise which randomly opens links 
between connected sites of the hypercube. The effects of 
this broken-link noise model have been previously consid- 
ered for a QW on a line [l^ and a on plane [1^ • The 
broken-link noise is an example of "unitary noise" j20j 
which may affect the dynamics in a different way from 
that resulting from performing coin or position partial 
measurements with probability p. This type of noise is 
characterized by a sequence of uncorrelated unitary op- 
erations applied on an initial state, 

|*(t)) = i7tC/t-i...C/i|*(0)) (16) 

and no measurements are performed during the evolu- 
tion. Each operator Ui, for z = l,...t, is of the form 
of Eq. (21), with a modified shift operator, S", which ac- 
counts for the current state of the network, i.e. which 
links are broken at time step i. 

The generalization of the shift operator should preserve 
unitarity. We take as a starting point the shift operator 
defined in Eq. 

n-l 2"-l 

where ej, a binary n-component vector with \ej\ — 1, 
labels the n spatial directions in the hypercube. The 
simplest way to generalize S to include the possibility 
of broken links is to define site-dependent, n-component 
binary vectors 

e'j{x)^e'j{x®ej) (18) 

{Bj, if link at site x in direction j is closed, 
0, if link at site x in direction j is open. 

Then, the modified shift operator is simply 

Ji-l 2"-l 

^'^E E \^®e'^ix)){x\®\j){j\. (19) 

If at site X the link in direction j is broken, 
|a; © e^-)(a;| = \x){x (B ej\ = \x){x\ and no probability flux 
is transferred across the broken link. The modified shift 
operator S' is unitary for any number of broken links. 
The evolution proceeds as follows. At each time step, 
the topology of the hypercube is defined, opening each 
link with probability p and setting the vectors e'j (x) ac- 
cording to Eq.dTH]). Then S"- (/(g) C) is applied to |*(t)) to 
generate the state at time t + 1. Note that in this model 
the state of the network at time < + 1 is uncorrelated with 
the previous state at t. 

In the presence of decoherence, the stationary distri- 
bution in the n-hypercube is independent of the initial 
conditions. In Fig. [6] we show the evolution of the total 
variation distance to (a) the uniform distribution and (b) 
the stationary distribution of the coherent case, Eq. 
It is clear that the introduction of even weak decoherence 



causes the asymptotic distribution to become uniform in 
a characteristic time ~ 1/p. After this time, the total 
variation distance to the uniform distribution decays ac- 
cording to an inverse power law which is independent of 
the decoherence rate p. For weak decoherence rates, the 
average probability distribution remains close to that of 
the coherent case for a time of order p^^/'^, as shown in 
panel (b) of Fig. [S] 

It is also apparent from Fig. [S] that the case p — 0.1 
goes faster to the uniform distribution than the other 
cases shown, both for smaller and larger values of p. In 
fact, when the decoherence rate is reduced or increased 
around this critical value, the convergence rate to the 
uniform distribution becomes slower. The fact that a 
critical decoherence rate Pc minimizes the mixing time 
has been previously reported by Kendon and Tregenna 
in the context of a decoherent QW in the iV— cycle. 
The dependence of the mixing time of the decoherent 
QW on the hypercube on the decoherence rate or, equiv- 
alently, the probability p of broken links is shown in the 
left panel of Fig. [T) A minimum can be identified near 
Pc ~ 0.1, which corresponds to a decoherence rate that 
may provide a faster mixing time in the hypercube. The 
critical value appears to be independent, or, at most, 
weakly dependent, on the dimension of the hypercube. 
This conclusion is similar to the one obtained by Kendon 
and Tregenna [15j for the A^-cycle with repeated mea- 
surements, although we have employed here a different 
kind of decoherence on a very different network. 

In the right panel of Fig. [7] we show the dependence on 
the dimension of the average mixing time to the uniform 
distribution for different decoherence levels. The linear 
mixing time to the stationary distribution, Eq. (|13p . of 
the coherent case is also shown (curve with circles). The 
decoherent mixing times increase with dimension at a 
rate which is approximately ri^^^ , i.e. slightly faster than 
quadratic, so that the mixing times with broken links 
are larger than the coherent mixing time for all dimen- 
sions. Note that the data in Fig. [7] also confirms that the 
walk with p « 0.1 mixes faster than the walks with other 
decoherence rates. 



IV. DISCUSSION 

In this paper the mixing properties of a discrete-time 
quantum walk on the n-dimensional hypercube has been 
considered in detail. The effect of decoherence from bro- 
ken links — a specific kind of unitary noise which involves 
no measurements — on the mixing properties of the walk 
has also been investigated. 

The stationary distribution for the coherent quantum 
walk on the rt-hypercube driven by a Grover coin has 
been found analytically, for the particular case of a sym- 
metric initial condition. It is not the uniform distribu- 
tion. However, the stationary distribution in this case is 
such that all Hamming weights are roughly equiprobable. 
It has also been noted that this distribution is dependent 
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Figure 6: Evolution of the total variation distance of the average distribution to (a) the uniform distribution, U = 2^" and (b) 
the stationary distribution, ■k{x) from Eq. (|12|) . Several decoherence rates are shown, together with the coherent case {p — 0), 
for fixed dimension n = 8. 
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Figure 7: Left panel: Average mixing time to the uniform distribution in a decoherent hypercube as a function of the probability 
of broken links p. Right panel: the same quantity as a function of dimension n. The average mixing time to the stationary 
distribution for the coherent case is also shown for comparison (curve with circles). 



on the initial condition. 

According to our numerical simulations the mixing 
time Me on the rt-hypercube increases as 0{n/e). We 
have shown this fact to be consistent with a general re- 
sult of Aharonov et al. [Tlj which, when particularized for 
the n-dimensional hypercube, provides an upper bound 
of 0(n3/2/e). 

The total variation distance to the uniform distribution 
shows a local minimum for t — ^n, as reported by Moore 
and Russell but a similar behavior with respect to 
the stationary distribution has not been observed. We 
have also found that the instantaneous mixing time to the 
stationary distribution, when it exists, has a nonlinear 
dependence on the dimension of the hypercube. These 
results reconcile the discrete-time quantum walk with the 
continuous-time quantum walk on the hypercube. The 
average distributions of both walks on the n-hypercube 
fail to converge to the uniform distribution, but both are 
uniform or e-close to uniform at certain times {t = jn). 

Decoherence, even at low rates, causes the stationary 



distribution to be the uniform distribution. The decay of 
the total variational distance takes place after a charac- 
teristic time 1/p and follows an inverse power law inde- 
pendent of p. For the case of broken-links unitary noise, 
which involves no measurements, an optimal decay rate 
for which the mixing time is a minimum has been found. 
In the mixing time as a function of dimension, the same 
optimal decay rate has been found to provide a faster con- 
vergence in the decoherent case. A similar effect has been 
reported by Kendon and Tregenna [l^ in the A^-cycle 
with decoherence from frequent partial measurements. 



In future work we are interested in generalizing these 
results in two ways: (i) finding the stationary distribu- 
tions induced by arbitrary initial coins and (ii) consider- 
ing different coin operations, specially the Fourier coin. 
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